Cross-fertilization of QCD and statistical physics: 
high energy scattering, reaction-diffusion, selective evolution, spin glasses 

and their connections* 
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High energy scattering was recently shown to be similar to a reaction-difTusion process. The latter 
defines a wide universality class that also contains e.g. some specific population evolution models. 
The common point of all these models is that their respective dynamics are described by noisy 
traveling wave equations. This observation has led to a new understanding of QCD in the regime 
of high energies, and known universal results on reaction-diffusion models could be transposed to 
obtain quantitative properties of QCD amplitudes. Conversely, new general results for that kind 
of statistical models have also been derived. Furthermore, an intriguing relationship between noisy 
traveling wave equations and the theory of spin glasses was found. 



I. INTRODUCTION 

The behavior of hadronic cross sections at high energy 
has been the subject of intense theoretical and experi- 
mental investigations for several decades. These studies 
will become even more relevant with the advent of the 
LHC, where protons will collide at a center-of-mass en- 
ergy of 14 TeV. 

In the following, we will focus on the forward elastic 
amplitude A{Y, k) for the scattering of two hadronic ob- 
jects at total rapidity Y, one of the objects (called the 
probe) being characterized by a tunable momentum fc, 
that defines the scales in the plane transverse to the colli- 
sion axis. The other hadron will be called the target. The 
scattering occurs at a given impact parameter that, here, 
will be fixed. Experimentally, the corresponding process 
could be, for example, deep-inelastic electron-proton (or 
nucleus) scattering, in which case k would be of order of 
the virtuality of the exchanged photon. The total rapid- 
ity Y can be seen as the rapidity of the target in the rest 
frame of the probe: It is related to the squared center- 
of-mass energy s through Y = ln(s/fc^). Note that one 
could equally well work in transverse coordinates instead 
of momenta, which proves useful in many cases: In coor- 
dinate space, A{Y, r) is directly related to the probability 
that the probe of size r interact with the evolved target 
at the considered impact parameter, which in particular 
results in the unitarity bound ^ < 1. The relationship 
between the two representations of the amplitude reads 



A{Y, r) = 



fk 



(1) 



Although total cross sections are not computable from 
standard QCD methods since confinement effects neces- 
sarily enter, the evolution with rapidity of the scattering 
amplitudes for a localized probe at a given impact pa- 
rameter (understood in Eq. ([1])) may be obtained starting 



from perturbative calculations that lead e.g. to the Bal- 
itsky hierarchy of equations 1] . In the latter framework, 
the evolution of A reads 



9ayA(y,r) 



2'K z^{r — z)' 



:iA{Y,z) + A{Y,r-z) 



A{Y,r)-{f{Y,z)f{Y,r-z))), (2) 



where a = — — - . One had to introduce T which is meant 
to be the scattering amplitude off one particular config- 
uration of the target. It is its average over all quantum 
fluctuations of the latter that corresponds to the physi- 
cally measurable amplitude: A — (T). 

Eq. (21) is not closed: One needs a further evolution 
equation for (TT). Eq. ([2]) is only the first equation 
of an infinite hierarchy which involves increasingly com- 
plex mathematical objects. Solving that system is a 
formidable task, and not surprisingly, no one has found 
a solution from a direct study of Eq. ([2]) . Another prob- 
lem with the hierarchy is that it is not clear whether its 
complete form for arbitrary targets is known. 

A significantly simpler equation may be obtained by 
arbitrarily factorizing the correlator (TT) into the prod- 
uct AA^ within a kind of mean field approximation that 
neglects fluctuations in the target. This was first pro- 
posed by Balitsky [l[ and subsequently re-derived by 
Kovchegov 0] in the particular physical context of deep- 
inelastic scattering off an infinitely large nucleus. An ele- 
gant representation of the resulting equation is obtained 
in momentum space by using Eq. ([1]), namely 



d^YA^x{-dink-)A-A^ 



(3) 
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where xi~(^\nk^) is the expression of the integral kernel 
that appears in Eq. Q in momentum space. The main 

Sropcrties of the solutions to Eq. Q are now known 
, JS, 6], but the effects of the fluctuations neglected in 
going to the Balitsky-Kovchegov (BK) equation ([3]) and 
even their very nature had not been appreciated until 
quite recently. 

One may guess that one difficulty with Eqs. ^ and ^ 
is that they both contain nonlinearities. The latter are 
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meant to encode parton saturation effects [?*] . One could 
think of neglecting them (by dropping the —A^ term 
in Eq. ([3|) for example), and indeed, that used to be 
the usual approximation until a few years ago. The 
resulting equation is named after Balitsky, Fadin, Ku- 
raev and Lipatov (BFKL) [§| {x{—dink^) that appears in 
Eq. jSl) is called the BFKL kernel). But in 1999, Golec- 
Biernat and WiisthofF showed 9] that the nonlinearities 
in Eqs. ([2]), ([3]) may already play an important role in the 
most recent data for deep-inelastic scattering (see also 
Ref. [l3|)- Their work stirred a great theoretical inter- 
est for the full nonlinear equations. For a long time, it 
was the BK equation that was the focus of most of the 
theoretical and phenomenological works in high energy 
QCD, despite the relative arbitrariness of the approxi- 
mations made to get it. Only in 2004 did Mueller and 
Shoshi address the problem of trying to quantify effects 
beyond those described by the BK equation, using a very 
original approach JTj. Subsequently, a deep interpreta- 
tion of what they had done in relation with problems of 
statistical mechanics was found 0, [l^, , which paved 
the way for a new and fruitful understanding of high en- 
ergy QCD. 

The goal of these two lectures is to show that high 
energy scattering looks very much like some processes 
(called reaction-diffusion) that appear in other physical, 
biological, or chemical contexts. We will provide a very 
simple and transparent picture of high energy QCD, but 
accurate enough to enable one to get presumably exact 
analytical results for QCD. The first lecture (Sec.|TT| aims 
at establishing this correspondence. "Cross-fertilization" 
is the title of these lectures because while investigating 
this relationship, we were able to also derive some general 
results that apply to a wider context. The second lecture 
(Sec. imp reviews them. The most important of these re- 
sults may directly be transposed to QCD. Throughout, 
the focus will be put on the underlying physical ideas 
rather than on the technicalities (references to original 
work where details are worked out will be provided when- 
ever necessary). We would like the reader to appreciate 
that this business was made possible by coming back to 
the basic concepts of high energy scattering in the par- 
ton model in the light of apparently completely unrelated 
physics, rather than trying to address directly the very 
technical equations that have been established for QCD. 



II. HIGH ENERGY QCD AS A STATISTICAL 
PROCESS 

A. What is universality? 

Our first task is to understand that it may sometimes 
be fruitful to replace a complicated problem (like QCD 
evolution, see Eq. ([2])) by a much simpler one (reaction- 
diffusion in our case, as we will see in these lectures), and 
yet be able to get quantitative results from the solution to 
the latter. This is possible when both models, the simple 



and the complicated ones, belong to the same universality 
class. 

The concept of universality was introduced in the 60's 
by a bunch of renowned physicists, among whom Widom, 
KadanofF, Fisher and Wilson in the theory of critical phe- 
nomena. To illustrate it, let us consider a system of spins, 
that is a set of binary variables Si = ±1, on a two- 
dimensional lattice and interacting through the Hamil- 
tonian 

i/ = -^J.,^.5,. (4) 

id 

If Jij — J > when the sites i and j are nearest neigh- 
bors and otherwise, this is the Ising model with fer- 
romagnetic interactions. The partition function is given 

by 

Z=Y, e-«[{^'>l/'=-^ (5) 

{S^} 

at a given temperature T. 

At high temperature, all spin configurations have equal 
probability. At low temperature, only a few configura- 
tions have a significant probability, the ones that mini- 
mize the energy. The minimum energy configuration is 
the one in which all spins are aligned. Near the critical 
point, the average magnetization reads 

(S,) ^ (T, - rf, (6) 

and an exact calculation gives (3 = 1/8 for two- 
dimensional systems, as may be checked in standard sta- 
tistical mechanics textbooks. 

The Ising model may represent only a very idealized 
magnet. One may wonder why it is interesting to study 
extensively such an over-simplified model, which obvi- 
ously incorporates only a few very basic properties of 
magnetic materials and neglects many others. Well, it 
turns out that the exponent P (as well as a few other 
exponents) is completely independent of the microscopic 
details of the system, which means that it is the exact 
result also for realistic magnets. Of course, other quan- 
tities strongly depend on those details: This is the case 
for the critical temperature Tc for example. 

The universality class is defined only by very general 
properties of the system, like the space dimension, and 
the symmetries. All systems that share these few prop- 
erties are expected to bear the same critical exponents. 
So if one is able to solve the simplest model, one is likely 
to know important features of all representatives of its 
universality class. 

High energy QCD has little to do with the Ising model 
and with criticality. We are now going to introduce a 
class of systems which, as we will argue later on, are likely 
to share universal properties with high energy QCD. 
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B. Warm-up: Properties of a stochastic 
zero-dimensional model 



We address the time evolution of a system of n{t) par- 
ticles that obey the following rules: In each time interval 
dt, each single particle has a probability dt to split in 2 
particles, and a probability — l)dt/N to disappear. 
This yields: 



n{t + dt) 



' n{t) + l proba dt 

/ X , , (n(t)-l)dt 
n{t)-l proba ^ ^ ' ' 



N 



(7) 



n{t) proba 1 — dt- 



{n{t)-l)dt 
N ■ 



One may represent the evolution of the system in the 
form of a stochastic evolution equation:-'^ 



dn n(n — 1) 

— = n 

dt N 



n{n — 1) 



N 



(8) 



where v is & noise (i.e. a random function of time) of zero 
average and of variance satisfying {v{t)v{t')) = S{t — t'). 
This means that v typically varies by 1 unit when t is 
changed by 1. It is a simple exercise to derive Eq. ([5]): It 
is enough to compute the mean {n{t + dt)) and variance 
(n^(t + dt)) given n{t) from the stochastic rules ([7]). 

We see on Eq. ([S])^ that in a typical realization, the 
number of particles will start to grow exponentially (lin- 
ear term) until n{t) N, at which point a steady state 
is reached (nonlinear term), up to random fluctuations 
(noise term). 

Solving the model consists in exhibiting the statistics 
oin{t). In a first approach, one may simulate numerically 
different realizations of the model, and then perform the 
average over the realizations to get {n{t)) and maybe 
higher-order moments of n{t). This is shown in Fig. [1] 
One could also try to obtain an evolution equation for 
{n{t)) from Eq. ([5]) and solve it analytically. But the 
latter involves the correlator (n'^): 



d{n) 
dt 



{n{n — 1)) 
N ' 



(9) 



and thus one eventually needs to solve an infinite hi- 
erarchy of coupled equations. This set of equations is 
qualitatively similar to the Balitsky hierarchy ([2]). Only 
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FIG. 1: (From Ref. Q^). Ten realizations of the evolution ^ 
of the zero- dimensional toy model (dotted lines), together 
with the full solution for {ii{t)) (full line) and the solution 
of the mean-field equation (|10|) (dashed line). A*' is set to 
5000 in this example. 



through a mean field approximation, consisting in factor- 
izing all correlators and valid in the large- iV limit does 
Eq. ^ boil down to a closed equation 



d{n) 
dt 



N 



(10) 



This is the equivalent of the BK equation ([3]). Unfor- 
tunately, the solution to this much simpler equation is 
quite far from the complete solution to Eq. ([8]) or (l9|) 
(see Fig. [T]), and it is manifest that more sophisticated 
approaches to the full stochastic model are needed. 

To address the complete stochastic model (O, we can 
think of two possible ways. One may compute iteratively 
the successive orders in 1 /N by field-theoretical methods 
and perform a resummation: This method was pushed 
quite far in Ref. 16]. Another possible path is to notice 
that the noise is actually important only as long as the 
number of particles n{t) is small: For large n(t), the evo- 
lution is essentially deterministic. For small n(t) instead, 
it is the nonlinearity in Eq. ^ that is negligible, which 
makes the evolution tractable, even analytically in this 
simple zero-dimensional case. That was investigated in 
Ref. [H. 



^ Eq. ((Sjl (as well as all other stochastic equations that appear in 
these lectures) has to be interpreted in the Ito way. 

^ Note that the form ^ is not particularly useful in practice, 
since i/ does not have any special properties. Our purpose in 
writing ^ is to get a feeling about the different possible repre- 
sentations of a stochastic process. We may also comment that 
there is a way to write the evolution as an Ito equation in which 
u would have a Gaussian distribution, and this is more interest- 
ing from a technical point of view. We refer the reader to [l5j 
for the details. 



C. Universality class of traveling wave equations 

We now add a spatial dimension to the model. We 
consider n{t, x) particles at time t and on each site x of 
a one-dimensional lattice. In addition to the rules ([7]) 
that apply on each site, a given particle at x has some 
probability p to jump to the nearby sites x+Ax or x—Ax. 
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The evolution of one realization may then be written as^ 



dn{t, x) 
Jt 



: p[n{t, X + Ax) + X — Ax) — 2n(t, x)] 



n(t,x)in(t,x) — 1) I — ; r 

^n{t,x)~ ^ ' ^' ^ l + ^f^^^v, (11) 

where v is again a noise of zero mean that varies by typ- 
icahy 1 when x or i are changed by one unit. This equa- 
tion is very close to Eq. ([8]) except for a new diffusion 
term (inside the square brackets) that correlates nearby 
spatial sites. 

Both in Eq. ([5]) and in Eq. (fTT|) the noise term is of 
order y^. This is typical of a statistical noise associated 
to discrete systems of particles: Indeed, adding stochasti- 
cally n particles on the average to a system means adding 
a number of particles typically in the range n ± ^Jn in a 
given realization. Hence the noise term is directly related 
to discreteness. 

Let us follow a particular realization of the evolution 
of the system. Starting from a given initial condition 
(for example a bunch of particles localized around some 
given lattice site), the number of particles grows exponen- 
tially on each occupied site until it reaches N ^ at which 
point the growth stops because the positive growth term 
n is exactly compensated by the negative nonlinear term 
—r?IN. At this point, only fluctuations of order \J~N are 
left in the evolution. At the same time, diffusion takes 
place, allowing the particles to "escape" toward larger 
(or lower) values of x (there is a symmetric front travel- 
ing to the left if the initial condition is local). It is then 
clear that the solution will look like a noisy wave front, 
that connects a region (say to the right) where there is 
no particle to a region (to the left) where the equilibrium 
number of particles per site N is reached. This wave front 
"travels" toward larger values of x, and is therefore called 
a traveling wave. One step of the evolution is illustrated 
in Fig. H 

The model that we have outlined is a typical reaction- 
diffusion process. The name "reaction-diffusion" stems 
from chemistry: Such a process may be a (very crude) 
model for chemical reactions. Renormalizing the number 
of particles through the introduction of u = n/N, one 
may write the general structure of the evolution of such 
a system in the following form: 

dfU — xi~dx)u — nonlinear function of u 

branching diffusion compensates the growth of u near 1 



noise of orders / — . (12) 



encodes discreteness 



^ To keep the equation simple and transparent, we replaced the 
noise term by its dominant part for n A'^, up to a constant. In 
other words, u has a variance that is only approximately normal- 
ized to unity. This does not change the discussion qualitatively. 





(a) 



(b) 



FIG. 2; Sketch of one step of the evolution of a realization in 
the one-dimensional model (from (a) to (b); A*' = 12 here). 
The dark disks represent the particles that have been added or 
that have diffused. For the smallest values of x, some particles 
have also disappeared, so that the maximum particle number 
does not exceed A'^ on the average, up to fluctuations of order 
y/N authorized by the stochastic term in Eq. (|lip . 



x{—dx) is a kernel that encodes branching diffusion. It 
could be a differential operator such as 9^ -I- 1 (which 
would correspond to carefully performing the limit Ax — > 
in the above example) or, more generally, an integral 
operator (which is the case in QCD). 

The universality class defined by Eq. (fT^ is usu- 
ally named after Fisher and Kolmogorov-Petrovsky- 
Piscounov (F-KPP) who for the first time for- 

mulated mathematically such processes. For a given 
stochastic equation, there is no general theorem to de- 
cide whether it belongs to the same class as Eq. (fT2)) 
or not. It is rather a matter of guess from the physical 
properties of the underlying model. One knowns how- 
ever that the exact form of the nonlinearity in Eq. (fT^ 
is not important. It could be any reasonable power of u, 
for example. Also the precise form of the noise is not an 
issue and does not affect the properties of the solutions 
at large N: One could have a slightly modified form, for 
example ^/u{l — u)/N. The only important feature is 
that it should scale like \/u/N for small values of u. 

We are now going to give some basic properties of the 
solutions to F-KPP-like evolution equations . More 
details and a review of the latest developments will be 
provided in Sec. IIIII 

The traveling wave that develops at large times is char- 
acterized by its mean shape and by the statistics of its 
position Xt. The latter may be defined, for example, as 
the value of x for which u has a definite value uo, for 
instance uq — 1/2. It is known since the seminal papers 
of Brunei and Derrida fl6| that for large times and large 
N, Xt has the following mean and variance: 



70 



2\ii^ N 



t 



(13) 



In'' TV 



The average is taken over many realizations of the evo- 
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FIG. 3: (From Ref. [T!]). 100 realizations of the evolution of 
a stochastic traveling wave model, at two different times. The 
averages over the realizations (that corresponds to the phys- 
ical amplitude in QCD, see the dictionary (|18p ) are shown 
in full line. Inset: the average (u) over the different realiza- 
tions, shifted in time and superimposed in order to show the 
mismatch in slopes. 

lution p2p . and xil) is the characteristic function of the 
diffusion kernel (i.e. the eigenvalue of x{~9x) that cor- 
responds to the eigenfunction e^''^). 70 solves x(7o) = 
7ox'(7o). The shape of each realization is, up to some 
noise 

M(t,a;) = e-'^o^^-^*). (14) 

A numerical evolution of a stochastic model of the F-KPP 
type is displayed in Fig. [S] One sees the universality in 
shape of each individual realization as well as the disper- 
sion in front positions between different realizations, that 
indeed grows like ^/^ (see Eq. We will come back 

to the derivation of these results in Sec. IIIII 



D. Evolution of QCD amplitudes 

With the intuition gained by studying simple particle 
models, we are now in position to argue that high energy 
scattering naturally belongs to the universality class of 
the stochastic F-KPP equation. 

We consider the scattering of two hadrons made of 
quark-antiquark pairs of respective sizes 1/k and l/fcg , to 
simplify the discussion (see Fig. [4]). We assume that these 
hadrons are small enough (much smaller than 1/Aqcd) 
so that perturbation theory applies. 

If there is few energy available for the scattering (the 
hadrons interact almost at rest), then typically they are 
in their valence configuration. Their interaction consists 
in exchanging one gluon when they are of similar size and 
at the same impact parameter: In this case, the forward 
elastic amplitude reads T ~ a^. If the hadrons are of 
very different sizes or if their impact parameters do not 



match, then they do not "see" each other and T = 0. 
Now boost one of the hadrons, called the target: This 
corresponds to increasing the center-of-mass energy of 
the reaction. Due to the resulting opening of phase space, 
a typical configuration of the target is no longer the bare 
valence quarks, but the latter together with a number 
n(y, k) of gluons of transverse momentum fc, whose mean 
depend on the rapidity. Correspondingly, the scattering 
amplitude off that particular configuration of the target 
reads 

T{Y,k)^aln{Y,k), (15) 

where k is of the order of the inverse size of the hadronic 
probe. This relationship is not exact, but may be derived 
rigorously as a limit of QCD scattering amplitudes.^ 

In a given realization of the Fock space of the target, 
the number of gluons typically doubles when aY is in- 
creased by one unit (each gluon has a probability to split 
in 2 gluons). The transverse momenta of the new gluons 
are close to the ones of their parents, up to some diffu- 
sion. This branching diffusion process is encoded in the 
BFKL evolution kernel x defined in the Introduction. As 
we are following one given realization, that evolves ran- 
domly, the equation contains a stochastic term of order 
-y/n, i.e. asVr from Eq. (fT5|) . At this stage, we may 
write the evolution of the scattering amplitude off one 
particular realization as 

davT ^x{-dink-)T + a,V^,^ (16) 

where is a noise of zero-mean whose variations are of 
order one when Infc^ and aY are varied by one unit. 

Because T is related to an interaction probability, 
T < 1 in appropriate normalizations and so the num- 
ber of gluons of given momentum (related to T) cannot 
grow exponentially forever. Hence a negative term that 
becomes important when T ~ 1 has to enter Eq. (|16p . 
One may write [H, [H, ^ 

daYT = x{-dink-}T-T^ + asV2Tiy. (17) 

The physical amplitude A is obtained from T by taking 
the average over all possible Fock realizations of the tar- 
get, A = (T). The mean-field approximation (T^) = A"^ 
of this equation precisely coincides with the BK equa- 
tion p. 

Hence we are led to an equation in the same universal- 
ity class as Eq. (fT^ . with the following dictionary (com- 



* Actually, it is useful to view the interaction in the framework of 
the color dipole model 20], in which gluons are assimilated to 
zero-size quark-antiquark pairs whose components combine with 
each other in the form of colorless dipoles, that evolve indepen- 
dently when rapidity is increased, n is in fact the number of such 
dipoles. 
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FIG. 4: One realization of the evolution of an onium-onium scattering amplitude. At zero rapidity (leftmost plots), the 
amplitude consists in a bare gluon exchange between the probe (of size 1/fc; upper bubble) and the target (of size l/fco; lower 
bubble): T = ou. li the momentum scales of the two objects match, T = Q otherwise. When rapidity is increased (going from 
the left to the right), each gluon in the Fock state of the target may split in two gluons of similar transverse momenta, and 
this process drives an exponential increase of the amplitude for k around fco, up to a diffusion. Eventually (rightmost plots), 
this growth has to slow down when T ~ 1 in order to comply with the unitarity constraints. Beyond this point, nonlinear 
saturation effects enter and a traveling wave forms in T. 



pare Eq. ^ and Eq. p7)) ): 

u ~ T, t = aY, a; = In fc^ , 

N ^ x{-dx) = BFKL kernel. ^'^^^ 

Let us comment on the physical implications of hav- 
ing a stochastic equation like (jl7p describing the dynam- 
ics. The BK equation ^ is known to admit traveling 
wave solutions [5] , whose position Xt (which corresponds 
in QCD to the so-called saturation scale Qs{Y) through 
Xt = In (5g(F)) is a deterministic function of the rapidity. 
This property is related to a phenomenon observed in the 
deep-inelastic scattering data and called "geometric scal- 
ing" [2l|: The physical amplitude A{Y,k) only depends 
on the combined variable k'^/Ql{Y), and not on k and Y 
separately. By contrast, each realization of an evolution 
given by Eq. ([T7]) is a noisy traveling wave, which has a 
random position (i.e. saturation scale) whose statistics 
follow Eq. (use the dictionary (HH])). The conse- 

quence is that A = (T) has a peculiar scaling behavior, 
different from geometric scaling [l^, : 



In fc2 - aY 



A{Y, k)^A 



( xilo) 
V 70 



7ox"(7o) 
21n^(l/Q2) 



^aY/\n\\lal) 



, (19) 



which is easy to derive from Eqs. ([13]) , (HH) . 



Mueller and Shoshi were the first authors who noticed 
that geometric scaling was broken beyond the BK equa- 
tion. However, the statistical interpretation given here 
was crucial to arrive at the correct scaling form ITOl) (see 
Ref. 03]). 

Let us summarize the assumptions that have lead to 
Eqs. pT|) and p^ . The linear terms therein represent 
the perturbative splitting of the partons when rapidity 
increases, and may be obtained exactly e.g. in the frame- 
work of the color dipole model [lO] ■ We have not written 
down explicitely the distribution of the noise v (we just 
gave its typical variations) , but the solution ([19]) is robust 
with respect to its detailed form. The nonlinear term is 
supposed to encode parton saturation. So far, we cannot 
provide an exact expression for it because a QCD real- 
ization of saturation has not been derived. However, the 
solution (fTO)) is also universal with respect to the form of 
the nonlinearities. The only really important condition 
for universality arguments to apply is that there should 
indeed be saturation, i.e. the gluon density should not 
be able to grow forever. Arguments for this feature were 
given a long time ago So if the latter is true, be- 
cause of the universality of some properties of solutions 
to Eq. p7|) . then we believe that Eq. p9)) represents the 
exact asymptotics of QCD. 

Last, we would like to draw the attention of the reader 
to the fact that we do not think that these results are 
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very relevant phenomenologically yet, since the asymp- 
totics show up only for In 1 /a^ ^ 1 (which means that 
as has to be unrealistically small). In our opinion, one 
may conduct sound phenomenological studies only once 
subleading effects have also been understood, beyond 
Eq. (I19p . However, numerical calculations may already 
provide some clues on the behavior of QCD amplitudes 
for more realistic values of the parameters, as explained 
in Ref. Ei. 



III. NEW GENERAL RESULTS ON 
TRAVELING WAVE EQUATIONS AND THEIR 
APPLICATIONS TO QCD 

In the previous section, we have argued that the evo- 
lution of QCD amplitudes at high energy is governed 
by a nonlinear stochastic evolution equation in the uni- 
versality class of the stochastic Fisher and Kolmogorov- 
Petrovsky-Piscounov (F-KPP) equation. The latter lit- 
erally reads 



dtu = d^u + u — u + \j ^"(1 ~ ")'^ 



(20) 



where is a normal Gaussian noise, uncorrelated both in 
t and in x. The deterministic part of this equation (ob- 
tained by setting ^ = 0) was first written in 1937 in the 
context of studies of the spread of genes (or diseases) in 
a population [l3|. Only in 1983 some properties of its so- 
lutions were rigorously derived p^ . In particular, it was 
understood that the large-time solutions were traveling 
waves. The dominant effect of the noise term in Eq. ([^D]) 
on the shape and velocity of the wave front was worked 
out by Brunet and Derrida in 1997. Note that the scal- 
ing form of the dispersion in the front position had only 
been measured numerically. We have recently been able 
to make progress on the statistics of the position of the 
front 

In practice, we will discuss the more general equa- 
tion (fT2l) . of which Eq. ((20)) is just a particular case. 



A. More on the propagation of noisy traveling 

waves 

1. Solving the deterministic F-KPP equation 

First, we address the deterministic version of Eq. (PO)) . 
namely 



dtu = x{~dx)u - 



(21) 



The choice x{~dx) = 9^ + 1 would correspond to Eq. (^0]) 
exactly, but we keep the branching diffusion terms in a 
more general form in order to be able to easily transpose 
the results to other contexts, such as the QCD evolution 
equations. 



In the case of the deterministic F-KPP equation, there 
are two fixed points, u = and u = 1 which are re- 
spectively unstable and stable under small perturbations. 
Starting from a localized initial condition, the branching 
diffusion encoded in x leads to a local exponential growth 
of u and a spread in x. When u becomes close to 1, then 
the nonlinear term tames the growth so that u does not 
get bigger than 1 . This mechanism leads to a wave front 
that invades larger values of x when time flows, being 
pulled by its tail. 

One may compute the asymptotic velocity of the front 
in a very simple way. Indeed, because the wave front is 
"pulled" by its low density tail, where the nonlinearity is 
negligible, it is enough to solve the linear part of Eq. (PT|) 
to understand the properties of the wave front. 



Consider a front decaying like 



7 a given wave number and v{^) its velocity. Then it is 
clear from Eq. pTjl that 



7 



(22) 



The general solution is a linear superposition of 

u{t,x) = / d7/(7)e-''(^~^('')*) 



'7: 



(23) 



where /(7) is a representation of the initial condition. 

Let us concentrate on the large-time behavior of the 
solution. We are following the wave around a fixed value 
of u (since the nonlinearity forces u to range between 
and 1), so the large-t limit has to be performed in 
the frame of the wave defined by the change of variable 
xwF = X — Voot, where Voc is the asymptotic velocity of 
the wave. After replacement of x in Eq. ([23]) . the saddle 
point condition reads Voo — x'(7)- Matching this expres- 
sion with Eq. (|22|) , we find that the wave number 70 that 
dominates asymptotically solves x(7o)/7o = x'(7o), i-e. 
7o minimizes 1^(7). Thus the asymptotic form of the wave 
front reads 



u(t c», a;) - e"''°("=-''('''')*) for x > v(7o)i, 

with f'(7o) — 0. 



(24) 



(The complete discussion may be found in Ref. 5] ) . This 
is the large time solution, but the front shape and veloc- 
ity are also known at subasymptotic times. Let Xi be 
the position of the front, defined e.g. by u(t, Xt) = mq, 
where uq is a constant between and 1. Then if the ini- 
tial condition is localized enough, the wave front has the 
shape [III 

u(t, x)^ix- Xt) exp ^^Mx-x.) (25) 

and the front velocity reads 

dXt x(7o) 3 



Vt 



dt 70 27ot 



(26) 
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The asymptotic velocity Voo (first term) is reached quite 
slowly. One more term is known in this large-t expansion 
0, but we will not discuss it here. We may comment 
that the velocity corresponds to the velocity of a 
front that has reached its asymptotic shape ((24|) only in 
the range 



(27) 



This is indeed the range in which the Gaussian factor in 
Eq. (P5|) is close to a constant. 



Hacking the deterministic equation to simulate 
discreteness 



Now we go back to the original stochastic equa- 
tions or (pn|) . Generally speaking, what we missed 
when we neglected the noise is mainly the discreteness of 
u. Think of a particle model on a lattice: u can take the 
values 2/7V, 3/N. . .but not fractions of these num- 
bers, a feature which of course was completely neglected 
in Eq. (HI]), as can be seen on the solution ([M)) . For 
u ^ 1/-/V, this might not be a problem, but for u ^ 1/N, 
neglecting discreteness may prove very wrong, especially 
since the propagation of the wave front is very sensitive 
to its tail. 

Brunet and Derrida came up with the idea that this 
discreteness could be incorporated back into the deter- 
ministic equation (|2ip by simply adding a cutoff^ that 
forbids the values of u between and [l^. They 
solved an equation of the type^ 



(28) 



and claimed that the solution matched the first order of 
the full equation in a small-noise (large- iV) expansion. 

We may estimate the large-i velocity Vbd of the trav- 
eling wave solution of Eq. . Starting from some ini- 
tial condition, the front evolves toward the asymptotic 
shape e"'''"'^^"'^'^ and its velocity increases according to 



^ This cutoff coincides with the one argued by Mueller and Shoshi, 
but the interpretation as being a consequence of discreteness of 
the number of partons had not been appreciated in Ref. Illl . 
Historically, we understood the relationship between QCD and 
statistical physics [T3l by looking for an interpretation of the 
Mueller-Shoshi results, in connection with the work of Brunet 
and Derrida [l9| . 

® Strictly speaking, the ©-function in Eq. (28) should not be ap- 
plied to the diffusion term. In the case of the F-KPP equation, 
one would rather write 



dtu ■■ 



although there would still be problems with the very mathemat- 
ical definition of such an equation. Eq. (28) would be correct 
literally if t were discrete with steps At of order 1. 



Eq. (1261) . However, as soon as the asymptotic front ex- 
tends down to values of u ~ 1/-/V, that is ioi x ^ Xt + L, 
where 



InA^ 
70 



(29) 



this shape cannot extend any further because of the cutoff 
in Eq. ((28)) . Then also the front velocity cannot grow 
any longer. According to Eq. ([27]) . this happens at time 
t ^ / (2x"(7o)), time at which the velocity reads T^bd = 
Vryo — 3cx"/(7o-^^^) according to Eq. (^5]) . c is a factor 
of order 1 whose determination needs a more accurate 
calculation [l9j. The complete result reads 



BD 



T^'^x"{la) ^ x(7o) 7i'^7ox"(7o) 
270-^^ 70 



2\iY^ N 



(30) 



The philosophy behind this procedure is quite transpar- 
ent. As soon as there are a few particles on a spatial 
site, parton number fluctuations are negligible and the 
subsequent evolution in that bin is deterministic. (One 
can get convinced of this statement already by looking 
at realizations of the evolution in the zero-dimensional 
model, see Fig. [T|). On the other hand, in the bins in 
which the deterministic evolution predicts a number of 
particles less than 1 (m < 1/N), one just sets m to in 
order to simulate discreteness. This is the effect of the 
cutoff. 



3. Incorporating stochastic effects 

The Brunet-Derrida cutoff procedure led to the fol- 
lowing result: The front propagates at a velocity Vbd 
lower than the velocity predicted by the mean-field equa- 
tion (HH), and its shape is the decreasing exponential 
(Eq. (211)), down to the position Xtip{t) = Xf° + lniV/7o 
{Xf^ = VbdOi which it is sharply cut off. 

But it may happen that a few extra particles are sent 
stochastically ahead of the sharp tip of the front (See 
Fig. [5]) . Their evolution would pull the front forward. To 
model this effect, we assume that the probability per unit 
time that there be a particle sent at a distance 6 ahead 
of the tip simply continues the asymptotic shape of the 
front, that is to say 



p{S) = Cie 



-70(5 



(31) 



where Ci is a constant. Heuristic arguments to support 
this assumption are presented in Ref. [22] (appendix A). 
Note that while the exponential shape is quite natural 
since it is the continuation of the deterministic solu- 
tion ((M)) . the fact that Ci need to be strictly constant 
(and cannot be a slowly varying function of S) is a priori 
more difficult to argue. 

Now once a particle has been produced at position 
a^tip + <5, say at time to, it starts to multiply (see Fig. [S]) 
and it eventually develops its own front (after a time of 
the order of L^), that will add up to the deterministic 
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Mean field evolution 
with cutoff 




the evolution of the position of the front Xt 



extra particle 



FIG. 5: Evolution of the front with a forward fluctuation. At 
time to, the primary front extends over a size L (Eq. (|29|l 'l and 
is a solution of the Brunet-Derrida cutoff equation (|28[) . An 
extra particle has been stochastically generated at a distance 
6 with respect to the tip of the primary front. At a later time, 
the latter grows deterministically into a secondary front, that 
will add up to the primary one. 



primary front made of the evolution of the bulk of the 
particles. 

Let us estimate the shift in the position of the front 
induced by these extra forward particles. Between the 
times to and t = to + L^, the velocity of the secondary 



front is given by Eq. ([26]). Hence its position 
relaxation, will be given by 



(2) 



xP^X^^ + 5+ f dt'vt, 

Jtn 



X 



BD 



to 

270 



after 



(32) 



where Xf^ = Vbd^- Eq. ([5^ holds up to a constant. 
We have used Eq. pS)) to express Vf-to- The observed 
front will eventually result in the sum of the primary 
and secondary fronts, after relaxation of the latter. Its 
position will be Xf^ supplemented by a shift R{S) that 
may be computed by writing the resulting front shape as 
the sum of the primary and secondary fronts: 



^^^^-,o(.-X--5+^lnL^)^ (33) 



where C2 is an undetermined constant. From Eq. (j33|) we 
can compute the shift due to the relaxation of a forward 
fluctuation R{S): 



1 / e'^°^ 
70 V 



(34) 



The probability distribution ([31]) and the front shift 
due to a forward fluctuation define an effective theory for 



Xt+dt — 



Xt + VBodt proba 1 — p{6)d6dt 
Xt + VBodt + R{S) proba p{6)dSdt. 



(35) 



From these rules, we may compute all cumulants of Xt: 



V-Vbb= J ddp{d)R{S) = 
[n-th ctimulant] 



C1C2 31nL 
70 70-^^ 

d5p{6)[RiS)r 

C1C2 n\({n) 



(36) 



70 76 



L3 



We see that the statistics of the position of the front 
still depend on the product C1C2 of the undetermined 
constants Ci and C2. We need a further assumption to 
fix its value. 

We go back to the expression for the correction to the 
mean- field front velocity, given in Eq. ([36]) . From the 
expressions of i?(^) (Eq. and of p((5) (Eq. dST])), we 

see that the integrand defining V — Vbd is almost a con- 
stant function of ^ for S < Sq = SIhL/jq, and is decaying 
exponentially for S > Sq. Furthermore, R{6o) is of order 
1, which means that when a fluctuation is emitted at a 
distance 6 ^ 60 ahead of the tip of the front, it evolves 
into a front that matches in position the deterministic 
primary front. We also notice that when a fluctuation 
has 6 < Sq, its evolution is completely linear until it is 
incorporated to the primary front, whereas fluctuations 
with 5 > So evolve nonlinearly but at the same time have 
a very suppressed probability. We are thus led to the nat- 
ural conjecture that the average front velocity is given by 
Eq. (ISni), with the replacement L lniV/70 + ^o- The 
large- iV expansion of the new expression of the velocity 
yields a correction of the order of In In N/ In^ N to the 
Brunct-Dcrrida result, more precisely 



V 



Xilo) 
7o 



7r^7oX"(7o) , 22//^ .31nlniV 

7oX (70) — riTTT- (37) 



2\n^ N 



70 In^ N 



Eqs. ((36)) and (|37|) match for the choice C1C2 = 
7r^x"(7o)- From this determination of C1C2, we also get 
the full expression of the cumulants: 



[n-th cumulant] 
t 



= ^^o'x"(7o) 



n!C(n) 
7j In^ N ' 



(38) 



We note that all cumulants are of order unity for t ^ 
In"^ N, which is the sign that the distribution of the front 
position is far from a trivial Gaussian, which makes it 
particularly interesting. On the other hand, they are 
proportional to k = t/\v? N, which is the sign that the 
position of the front is the result of the sum of k indepen- 
dent random variables, and as such, becomes Gaussian 
when K is very large. 

These results rely on a number of conjectures that no- 
one has been able to prove so far. However, we performed 
very precise numerical checks on specific models, and we 
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found a perfect matching (see Ref. [2^). So we are rea- 
sonably confident that our expressions are the correct 
ones. Providing a formal proof would be an interesting 
challenge for a mathematician. 

The analytical results presented here, namely Eqs. (|37p 
and ([38|) . may directly be applied to QCD using the dic- 
tionary (|18p . This would help to improve the analytic 
form of the amphtude (fTT]) . 



B. Traveling waves in evolution models 

Let us now consider a model of population evolution. 
Each individual is characterized by a single real number x 
that measures its adequacy to the environment. At time 
t, there are n{t, x) individuals in the population with a 
given adequacy x. To go from time t to t -f 1, we give 
the following rule: Each individual dies after having given 
birth to 2 offspring, that have respective adequacies x+ei 
and X + £2, where ei,2 are random variables distributed 
according to a sufficiently local probability distribution 
tp. If the total population exceeds a given integer iV, then 
we get rid of the individuals with the lowest values of x in 
order to keep the population size constant and equal to 
N . This model could represent for example the evolution 
of a population of bacterias under asexual reproduction, 
in a medium where resources are limited, which enforces a 
selection of the "best" individuals. Let us define u{t^ x) 
to be the fraction of population that has its adequacy 
variable larger than x. It is not difficult to realize that u 
has the shape of a wave front, connecting a region where 
u — 1 (for small values of x), and a region where u = 
(for larger values of x), that moves toward positive x 
when time elapses, see Fig. [HI 

One may write a stochastic evolution equation for u: 

u{t + 1, a;) = min (1,2/ de'ip{e)u{t, x — e)) + 



J de4}{e)u{t,x-e)iy{t + l,x) (39) 

where as usually, v is a, noise of zero mean and variance 
of order 1 in the region where u is small. At first sight, 
it is far from obvious that this equation should be in the 
universality class of the F-KPP equation. First, this is 
a finite difference equation in t rather than a differential 
equation. Second, the noise term is non-Gaussian and 
furthermore, a closer look shows that v is strongly cor- 
related spatially. The nonlinearity is also of completely 
different nature in this model. However, the fact that 
a noisy front is formed for u and especially the physical 
mechanism of evolution (branching diffusion-f saturation) 
points toward the F-KPP universality class. If this is 
true, we can apply the results obtained in the previous 
subsection to study how the mean adequacy x grows with 
time (this would be the average position of the front). It 
is easy to see that for this model, the characteristic func- 



X 




FIG. 6: Evolution of the population over a few generations. 
In this model, each individual of a given generation t has ex- 
actly two offspring, whose adequacy to environment differs 
from his by the random variables ei and £2. Once the pop- 
ulation size has reached (A^ = 8 in this sketch), it is kept 
constant by allowing only the "best" individuals to have a 
descendance. The fraction u of individuals having their ad- 
equacy variable larger than x is plotted (bottom plot) for 3 
different times. One sees that u exhibits the characteristics 
of a traveling wave. 

tion would be 

X(7) = 111 (2 / rfeV'(£)e^") . (40) 

We checked that a numerical simulation of this model in- 
deed matches the analytical predictions, confirming that 
the model (|39p is indeed in the universality class defined 
by the F-KPP evolution. 

For that kind of evolution models, one may ask a fur- 
ther question. At a given time t, one may pick a given 
number of individuals, say 2 or 3, and count the number 
of generations one has to go back in the past to find their 
most recent common ancestor. We denote these numbers 
by T2, T^. ■ ■ This question may be relevant to studies of 
the genetic diversity of a population. 

It is not difficult to guess the order of magnitude of 
Ti. We saw in the previous subsection that Xt looked 
like the sum of k = t/ln^ N independent random vari- 
ables. One can understand this fact very precisely on this 
model. Indeed, from time to time, the evolution gener- 
ates an individual whose adequacy to the environment is 
much larger than the one of any other individual. His 
offspring will partly inherit his adequacy, according to 
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FIG. 7: Evolution of the population over a few generations. Each individual has two offspring in the next generation, but only 
the A*' 8 in this illustration) with highest values of x are allowed to have a descendance. Sometimes, an individual is so 
advanced (it has to be ahead by at least 5 = 5o = In In N) that his descendance will take over the whole population. This 
happens, on the average, once in In'' generations. An estimator of this time scale may be obtained by picking 2, 3 or more 
individuals at random and looking for their common ancestor. 



the stochastic evolution rules and thus will still be in ad- 
vance with respect to the bulk of the population. After 
some generations, only his descendants may survive, so 
that he will be the common ancestor of the whole popu- 
lation. For an individual to have a significant probability 
to have his descendance replacing the whole population 
at some later time, the position of the secondary front 
developed by this individual has to be larger than the 
position of the primary front, a condition which reads 
5 > lnL^/70. According to Eq. ((3T|) . this happens once 
in = (lnA^/7o)"^ generations. Hence 

T, ~ In^ N. (41) 

It turns out that through a slightly more elaborate anal- 
ysis, we may also fully compute the coefficients of the Ti, 
simply from the assumptions of the previous subsection 
on the mechanism of the propagation of the front. We 
found m 

^ 5 (74^ ^ 25 
(T2) 4' {T2) 18' ^ ' 

These numbers characterize the statistics of the genealog- 
ical trees that can be drawn in this kind of models. 

It is not clear if results such as (|41]l and (j42|) may be 
applied to particle physics observables. The latter may 
be more relevant to evolutionary biology. 



C. Spin glasses 

A spin glass [131 is a system of spins with random inter- 
actions. The Hamiltonian is given by Eq. but the Jy 



FIG. 8: Two frustrated configurations of 4 spins, that have the 
same energy. The and "— " on the lattice links represent 
the values of the interaction Ji, . 



are now random numbers endowed with a given probabil- 
ity distribution. The simplest case is when Jy are binary 
variables ±1, but the standard choice is a Gaussian dis- 
tribution which leads to a model named after Edwards 
and Anderson. The Sherrington-Kirkpatrick model is the 
infinite-range version of the latter: The same distribution 
holds for any pair of spins and not only for nearest neigh- 
bors. 

This time, the minimum energy state at zero tempera- 
ture is not unique, contrarily to the Ising case for which 
the minimum energy configuration is reached when all 
spins are aligned. This is due to "frustration" (the fact 
that all interactions can never be satisfied simultane- 
ously), and is illustrated in Fig. [51 

To classify the configurations, one may define the over- 
lap between two spin configurations a and 6 as a function 
that reflects the fraction of spins that arc common to a 
and b. The usual definition is qab = J2i ^f^i- Speak- 
ing through the hat, the configurations that minimize 
the energy are organized ultrametrically with respect to 
a distance derived from the overlap. This means that 
the configurations may be represented as the leaves of a 
(random) tree. The statistics of this tree, first found by 
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Parisi [23] , turn out to be exactly the same as the ones of 
our genealogic trees in the previous subsection that are 
encoded in the ratios (j42| . 

It should be clear to the reader that there is no a priori 
reason for such a feature. The method for finding the 
statistics in Parisi's theory relies on the replica method 
(see Ref. [1^ for details) , and has nothing to do with the 
way we derived the result in our case. Whether this link 
is deep or accidental is an outstanding open question, for 
which we have no clue. 



IV. CONCLUSION 

When it was first suggested that high energy QCD has 
something to do with reaction-diffusion processes and 
that this analogy leads to new quantitative predictions 
[l^ . many experts were skeptical and reluctant to ac- 
cept this new point of view. Getting these ideas through 
has involved some passionate debates, to say the least. 
But at the time of the present lectures, that is two years 
later, this proposal has undoubiously triggered a burst of 
activity, that has spread in many different directions. 

The main advantage of viewing high energy scattering 
in the peculiar way that we have explained here is that it 
provides a simple picture endowed with a clear physical 
interpretation. This contrasts with the high technicity of 
evolution equations for QCD amplitudes that had been 
derived over the last 10 years such as the Balitsky hi- 
erarchy ([2|). Furthermore, this picture has a direct and 
intuitive connection with the basics of the parton model. 



This approach has paved the way for new results in QCD, 
obtained from tools known in statistical physics. 

Obviously, there is still room for improvements of our 
understanding of high energy scattering. Admittedly, our 
proposal relies on a few conjectures that will have to be 
proved. In our opinion, the most problematic one is the 
assumption that the number of gluons or of color dipoles 
saturates at l/a^ per unit of transverse phase space. Of 
course, this has been part of the "folklore" of high energy 
physics at least since the seminal work of Gribov, Levin, 
Ryskin in the 80's Q- However, no QCD realization of 
saturation has yet been found. A relevant question would 
be, for example, how parton recombination may be im- 
plemented in the framework of the color dipole model. 
Some interesting ideas have recently been put forward 
(see e.g. Ref. [26|), but a rigorous proof in QCD still has 
to be provided. 

In the second lecture, we have presented some recent 
advances in statistical physics to which we have been 
able to contribute. Here also, there is an intriguing 
correspondence between two seemingly different physical 
problems: stochastic fronts and the theory of spin 
glasses. Whether this is something deep or simply 
accidental still remains to be understood. 
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